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Abstract 

We present a method for cleanly extracting the CP phase 7 from the Dalitz plots 
of B+ K+tt+tt-, B° d K + 7i°7r-, B° d -> K®7r~^7t~ 7 B% -> K + K°K~, and 5^ -> 
K°K°K°. The B ^tttt and 5 -»■ decays are related by flavor SU(3) 

symmetry, but SU(3) breaking is taken into account. Most of the experimental 
measurements have already been made - what remains is a Dalitz-plot analysis of 
B d — > K°K°K° (or B d — > KsK$Ks ). We (very) roughly estimate the error on 
7 to be ~ 25%. This is somewhat larger than the error in two-body decays, but 
it would be the first clean measurement of 7 in three-body decays. Furthermore, 
at the super-S factory, it is possible that 7 could be measured more precisely in 
three-body decays than in two-body decays. 
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In the past, most of the theoretical work looking at clean methods for extract- 
ing weak-phase information in the B system focused on two-body decays. This is 
essentially because (i) final states such as ipKs, tt + tt~, etc. are CP eigenstates, and 
(ii) if there is a second decay amplitude, with a different weak phase, it has been 
possible to find methods to remove this "pollution," and cleanly get at the weak 
phases. On the other hand, in three-body B decays, final states such as Kstt + tt~ 
are not CP eigenstates - the value of its CP depends on whether the relative tt + tt~ 
angular momentum is even (CP +) or odd (CP — ). Furthermore, even if the CP of 
the final state were determined in some way, one still has the problem of removing 
the pollution due to additional decay amplitudes. For these reasons, it has generally 
been thought that it is not possible to obtain clean weak-phase information from 
three-body decays pQ. 

Recently, it was shown that this is not true. By doing a diagrammatic analysis of 
the three-body amplitudes, one can resolve these two problems [2]. First, a Dalitz- 
plot analysis can be used to experimentally separate the CP + and — components 
of the three-particle final state. Second, one can often remove the pollution of 
additional diagrams and cleanly measure the CP phases. In fact, in Ref. [3J, it was 
shown how to extract the weak phase 7 from B — > Ktttt decays. We briefly describe 
this method below. 

In B — > Ktttt decays, the isospin state of the tttt pair must be symmetric (anti- 
symmetric) if the relative angular momentum is even (odd). As we will see below, it 
is the symmetric case which is most interesting. Here there are six possible decays: 

B + K + TT + TT-, B + K + TT TT , B + -> K°TT + TT , B° d -> K + TT~ TT° , B° d -> K°TT + TT-, 

and B® — > K°tt tt°. The first step is to express the amplitudes for these processes in 
terms of diagrams. The diagrams are as in two-body B decays [I] : the color-favored 
and color-suppressed tree amplitudes T and C, the gluonic-penguin amplitudes P tc 
and P uc , and the color-favored and color-suppressed electroweak-penguin (EWP) 
amplitudes Pew an d Pew- (We neglect annihilation- and exchange- type diagrams.) 
Furthermore, for three-body decays, it is necessary to "pop" a quark pair from the 
vacuum. The diagrams are written with subscripts, indicating that the popped 
quark pair is between two (non-spectator) final-state quarks (subscript '1'), or be- 
tween two final-state quarks including the spectator (subscript '2'). (For B — > Ktttt 
decays, the popped quark pair is uu or dd. Under isospin, these amplitudes are 
equal.) 

In addition, some time ago it was shown that, under flavor SU(3) symmetry, 
there are relations between the EWP and tree diagrams in B — > Kit decays [51 E] • 
In Ref. [3J , it was shown that similar EWP-tree relations hold for B — > Ktttt decays. 
Taking C\/cz = Cg/cio for the Wilson coefficients (which holds to about 5%), these 
take the simple form 

PeWI = > PeW2 = K ^2 J 

PeWX = ' PeW2 = K ^2 ' (1) 
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where 



3 |A t | Cg + Cio , v 

K = ~o^b — T — ' 2 
2 |Ai SJ | ci + c 2 

with a« = f; 6 v ps . 

Now, the EWP-tree relations assume SU(3) symmetry (and the approximate 
ratio of Wilson coefficients). The expected error due to SU(3)-breaking effects is 
0(30%). However, the dominant diagram in b — > s decays is P/ c , so that EWPs 
and trees are subleading effects. Thus SU(3) breaking is subdominant - the net 
theoretical error due to the use of the EWP-tree relations is only 0(5%). This 
is consistent with the error estimates given in Ref. [5] (for EWP-tree relations in 
B ->■ Kir). 

In addition, there is an important caveat. Under SU(3), the final state in B — > 
Kirn involves three identical particles, so that the six permutations of these particles 
(the group S3) must be taken into account. That is, the three particles are in a totally 
symmetric state, a totally antisymmetric state, or one of four mixed states. However, 
the EWP-tree relations hold only for the totally symmetric state. Thus, the analysis 
must be carried out for this state. Now, the expressions for the B — > K(n7i) sym 
amplitudes in terms of diagrams hold even under full SU(3) symmetry [3J. It is 
therefore only necessary to produce observables for the totally symmetric states. 
This is doable, and below we present the details of how this is carried out. 

With the above EWP-tree relations, the six B — > K(im) S y m amplitudes can be 
written in terms of 5 effective diagrams (i.e. linear combinations of the diagrams) 
[3j. There are therefore 10 theoretical parameters in the amplitudes^: 5 magnitudes 
of effective diagrams, 4 relative (strong) phases, and 7. On the other hand, there 
are 11 experimental observables. Given that B + — > K°tt + tt° is not independent (its 
amplitude is proportional to that of B d — > K + tt°'k^), these are the branching ratios 
and direct CP asymmetries of B + -> K+tt+tt', B + -> K + n n , B° d -> K + 7i 7r~ , 
B d -> K°7i + 7i~, and B d -> K°n ir°, and the indirect CP asymmetry of B d -»■ 
K°7t + 7t~ (the indirect CP asymmetry of B d — > K 7t°tt° will essentially be impossible 
to measure). Since there are more observables than theoretical parameters, 7 can 
be extracted by doing a fitQ. 

The disadvantage of this method is that it involves the decays B + — > K + ir 7r° 
and B d — > K°tt 7i . With two 7r° mesons in the final state, both of these decays 
will be extremely difficult to measure. We are therefore motivated to see if the 

3 In fact, the expression for any indirect CP asymmetry contains another theoretical parameter 
- the phase of B^-B d mixing, j3. However, its value can be taken from the indirect CP asymmetry 
in B° d J/^Ks m- 

4 There is a complication in that the diagrams are momentum dependent, as are the observables. 
In obtaining the best-fit "values" of the diagrams, one will determine the momentum dependence of 
their magnitudes and relative strong phases. On the other hand, 7 is independent of the particles' 
momenta. Later in the paper, we detail how such a fit is done. 
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P — > Ktttt method can be modified, avoiding these two decays. As we show below, 
this can indeed be done - things can be considerably improved by using P — > KKK 
decays. The use of these decays is quite natural since they, like B — > Ktttt, are also 
b — > s transitions. 

First, consider B — > Ktttt decays with the tttt pair in a symmetric isospin state. 
We leave aside B + K + tt°tt°, B° d K tt tt and B + K n + n (since, as men- 
tioned above, its amplitude is not independent). The amplitudes of the remaining 
three processes are 

— * K + tt°tt ) syrn = Tje* 7 + C 2 e t ' y — P^v/2 ~ Pewi > 
V2A(B° d -+ K\ + K-) sym = -T[e^ - C[e^ - P' uc e^ + P[ c 



1 p' _l 2 P /C _ P'C 

i* ~^ r EWl ~^ r EWl ~^ r EW2 > 

V2A{B + ^ K + Tj + Tj-) sym = — T^e* 7 — C[e ij — P' uc e n + P' u 



tc 

4- - p' --p ,c 4- - p /C m 

^ ^ r EW\ <^ r EW\ ^ ^ jr EW2 ■ \°J 

In the above, P' = P{ + P^. (As P — )■ i^7T7r is a 6 — >■ s transition, the diagrams are 
written with primes.) Here we have explicitly written the weak-phase dependence 
(this includes 7 and the minus sign from V^ b V ts [P/ c and EWPs]), while the diagrams 
contain strong phases. 

Second, consider P — > KKK decays. For the case in which the final KK pair 
is in a symmetric isospin state, there are four such processes: P + — > K + K + K~ , 
P+ ->• K + K°K°, Bj K + K°K~, and P^ ->• K°K°K°. Here, P+ -4- K+K + K~ 
and P + — >■ K + K°K° are not independent - their amplitudes are proportional to 
those of B Q d ->- K + K°K- and P^ -4- K°K°K°, respectively. These are 

>/2A(Bjj ^ + ir ir-) s , m = -T^ 7 - C( iS e i7 - P^ 7 + P/ c 

, 2 p/ 1 p/ 1 2 p/C 1 p/C 

A(P° -> K°K°K°) sym = PLe^-PL (4) 



J p' _| p' _| p'C _| p'C 



where P' = Pg s + P[. In the above, certain diagrams are written with the subscript 
£ s.' This indicates that the popped quark pair is ss. When the diagram has no 
subscript s (the penguin or EWP diagrams), this means that the popped quark pair 
is uu or dd, but the virtual particle decays to ss. 

We now assume flavor SU(3) symmetry. This has two consequences. First, the 
amplitude with a popped ss quark pair is equal to that with a popped uu or dd. That 
is, we no longer need the subscript s on diagrams. This means that the diagrams 
in P — > KKK decays are the same as those in P — > Ktttt decays. Second, the 
EWP-tree relations of Eq. (JTJ hold. 
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Thus, under SU(3) the amplitudes of Eqs. and (HI) take the form 

2A(B d ^K + n 7r-) sym = T[e* + C 2 e* -k{T^ + C[) , 

V2A(B° d ^K% + n-) sym = -T^-C'^-P^ + Pi 

f\ 2 1 
+ K ( Lt[ + -C' - -C' 
V3 1 3 1 3 2 

V2A(B + -> K + n + n') sym = -T^-C[^-P' uc ^ + Pi 

/1 1 2 

+ k ( -t: - -a, + -a 

V3 1 3 1 3 2 
V2A(B° d -)■ K + K°K-) sym = -Tie* - C[e* - P' uc e* + P/ ( 



Note that this implies that /\ + 7r+7r~) S3/m = A{B Q d K+K°K-) sym . Fur- 

ther, we reiterate that the above expressions for the amplitudes hold also for the 
totally symmetric final state, to which the EWP-tree relations apply 
We now define the following five effective diagrams: 



r 

a 


= T[-Tl, 


T' 


= C 2 4" T2 1 


P' 

a 








H 






= k{C[-C 2 ) . 



3 1 3 2 

(6) 

The amplitudes can be written in terms of these five diagrams: 

2A(B° d -+K+n\-) sym = T a ^ + T^-C' a - K Tl, 

V2A(B° d ^K\+7r-) sym = -r a e* - I*e* + I* b , 

V2A(B+ -> K + n + n-) sym = -P' a e* + P' h - C' a 

V2A(B° d ^K+K°K-) sym = -P^ + Pl-C' a 

A(B° d ^K K K°) sym = P' a e* -T' h e* ~\c' a e* 

- Pl + K T a + KTl + C' a . (7) 

As with the B — > Ktttt method, five effective diagrams corresponds to 10 theo- 
retical parameters: 5 magnitudes of diagrams, 4 relative phases, and 7. But there 



a 1 

'a 1 
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are 11 (momentum-dependent) experimental observables: the decay rates and direct 
asymmetries for the four decays B d — > K+n ^, B d — > /^°7r + 7r~, B d — > K + K°K~ 
and B d — > K°K°K° (we ignore B + — > K + -k + -k^ since its amplitude is not inde- 
pendent), and the indirect asymmetries of B d — > K 7t + tt~, B d — > K + K°K~ and 
B d K°K°K°. With more observables than theoretical parameters, 7 can be 
extracted from a fit. 

We now present the details of how the fit is carried out. Consider the decay 
B —7- P1P2P3, in which the three pseudoscalar mesons Pi (i = 1-3) have momenta 
Pi. From these, we can construct the three Mandelstam variables: 

S12 = (pi +P2) 2 , S13 = (Pi +P3) 2 , s 23 = (P2 + P3) 2 ■ (8) 
These are not independent, but obey 

Sl2 + «13 + -§23 = m B + m i + m l + m i ■ (9) 

Experimentally, the Dalitz plot of this decay is measured. Its events are given in 
terms of two Mandelstam variables, say si 2 and S13. Now, the great advantage of a 
Dalitz-plot analysis is that it allows one to extract the full amplitude of the decay. 
We write 

M(B -> P!P 2 P 3 ) = J £c j e^F j (s 12 ,s 13 ) , (10) 

3 

where the sum is over all decay modes (resonant and non- resonant). Cj and 9j are 
the magnitude and phase of the j contribution, respectively, measured relative to 
one of the contributing channels. The distributions Fj, which depend on s 12 and s 13 , 
describe the dynamics of the individual decay amplitudes, and take different (known) 
forms for the various contributions. The key point is that a maximum likelihood fit 
over the entire Dalitz plot gives the best values of the Cj and 6j. Thus, the decay 
amplitude can be obtained, up to an overall normalization. This normalization is 
fixed by the constraint of the measured partial rate [7] : 

r = / l"^.^3. (ID 

With this, the decay amplitude M.(s 12 , S13) is known. 

As will be seen below, we rely heavily on -M(si2, S13). In particular, we use it 
to obtain the observables for the B — > P1P2P3 decay. As such, the errors on these 
observables come entirely from the uncertainty in (si 2 , S13). While, as noted 
above, it is possible to obtain the best-fit values of the Dalitz-plot variables Cj and 
9j, there are errors associated with these values. This is due to two sources. First, 
one has the statistical error in the experimental Dalitz plot. Second, there is a 
systematic uncertainty related to the choice of the Fj in Eq. (flOl) . In addition, there 
is a statistical error in the overall normalization [coming from Eq. fllip ] . All of these 
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must be carefully taken into account in order to obtain conservative errors on the 
Dalitz-plot variables. 

As noted earlier, the EWP-tree relations hold only for the totally symmetric 
SU(3) decay amplitude. But this can be found from the above: 

M fully sym = [M (s 12 , S 13 ) + M (s 13 , S 12 ) + M (s 12 , S 23 ) 

+ M{s 23 ,s 12 ) + M(s 23 ,s l3 ) + M{s 13 ,s 23 )} . (12) 

Using this, it is possible to compute the B — > P 1 P 2 P 3 observables. However, recall 
that the method involves a fit using the observables from several different decays 
{B° d K + n°it-, B° d -»■ K°ir + n-, B° d -»■ K + K°K~ and B° d K°K°K°). All 
observables must involve the same Mandelstam variables. On the other hand, the 
numbering of final-state particles is arbitrary, so that Si 2 for one decay might equal 
S13 for a different decay. All of this makes it somewhat confusing to ensure that 
observables in different decays have the same Mandelstam variables. For this reason, 
it is useful at this stage to change notation (but the physics is unchanged). In any 
decay there are three Mandelstam variables. We define s ++ , s + and s_ to be the 
largest, second-largest, and smallest of these, respectively. The identities of the 
particles which are associated with s ++ , s + and s_ are irrelevant (e.g. s ++ can 
correspond to Si 2 , S13 or s 23 ). This is consistent with the assumption of SU(3) and 
the fully symmetric decay amplitude. With these Mandelstam variables, we have 

M fully sym = [M (s ++ , 8+) + M (s+, S ++ ) + M (s++, S_) 

+ M(s„,s ++ ) + M(s.,s + ) +M(s+,s-)] . (13) 

Since s ++ , s + and s_ are not independent, this gives the fully symmetric amplitude 
as a function of two Mandelstam variables, say s ++ and s+. 

The observables are obtained as follows. First, one forms the totally symmetric 
SU(3) decay amplitudes as in Eq. ([13]) for each B — >■ PiP 2 P 3 decay (M.f u u y sym ) and 
its CP conjugate (M.f u u y S ym)- Second, using these, for specific values of s ++ and 
s + , one computes the partial rates: 

1 3- I -M fully sym(s++, S+)\ 2 



(2vr) 3 32m| 



r s++ , s+ - ( 27r )3 32m | \Mfvllv sym(s ++ , S + )\ 2 . (14) 

These allow the computation of the CP-averaged branching ratio and direct CP 
asymmetry: 



r( r - 

J- B V 



r — r 
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Third, for those decays in which the final state is accessible to both B d and B d 
mesons, one has an indirect (mixing-induced) CP asymmetry. It is given by 



e 



-2i/3 



M fully sym S+) 
M fully sym\£>-\ — h "+/ 



(16) 



As discussed earlier, in all cases, the error on the observables is found by propogating 
the errors on the Dalitz-plot variables. These include both statistical and systematic 
effects. 

Now, given that the method assumes flavor SU(3) symmetry, one would like to 
know how SU(3) breaking affects the analysis, and what is its size. Leaving aside the 
EWP-tree relations, in which SU(3)-breaking effects are subdominant, there are two 
areas where the breaking may be significant. First, under SU(3), the diagrams in 
B — > KKK and B — > Ktttt are the same. Since both decays are b — > s transitions, 
the difference between them is that B — > KKK decays have an ss quark pair in the 
final state, hadronizing to KK, while B — > Ktttt decays have uu or dd, hadronizing 
to 7T7T. This is essentially the same for each diagram. (The SU(3)-breaking effect 
associated with an ss pair being popped from the vacuum may not be exactly equal 
to that when ss is produced in the decay of a virtual particle, but the difference is 
small.) Thus, including SU(3) breaking, the amplitudes of Eq. ([7]) can be written 

2A(B° d -+ K+n\-) sym = T^ + T^-C f a -Kn, 

V2A(B° d ^K°7r + n-) sym = -T' a e^ - P' a e^ + P 6 ' , 

y/2A(B+ -+ K + n + TT-) sym = -P> a e*r + PZ-C' a , 

V2A(B° d ^K + K°K-) sym = (l + f sm )\-Pyy + Pl-C' a 



A(B° d ^K K K°) sym = (l + jW (3) ) 

- Pl + K T a + KTl + C' a ] , 117) 



P'e 47 - Tie 1 ' 1 - -C'e^ 



where fsu(3) is the SU(3)-breaking factor. Second, under SU(3), 7r's and K : s are 
identical particles, so that there is no difference between the Mandelstam variables 
for the processes B — > KKK and B — > Ktttt. There is therefore an SU(3)-breaking 
effect between the fully symmetric decay amplitudes for the two types of decay. 
However, it can be included in fsu{3)- 

The addition of fsu(3) brings the number of unknown theoretical parameters 
to 11. In principle, these can all be determined from a fit to the 11 experimental 
observables, albeit with discrete ambiguities. However, we can do better. Above it 
was noted that, in the limit of perfect SU(3), A(B + — > K + TT + TT~) sym = A(B d — > 
K + K°K~) sym . This means that fsu(3) can be determined by a comparison of these 
two decays. In particular, 



t + B(B d — >■ K + K°K~) 
B(B+ K+tt+tt-) 



^ = (i + w 

sym 
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In fact, this comparison can be performed now since the decays have been mea- 
sured: B + K+7T+7T- in Ref. [BJ, B%-> K + K°K~ in Ref. [9]. Now, since the EWP- 
tree relations [Eq. ([I])] have been used to derive the expressions for the amplitudes, 
Eq. ( 118)) holds only for totally symmetric states. Using the technique described 
above, one can obtain A(B + -y K + n + 7i') fuUy sym and A{B° d -)> K + K K')f uUy sym . 
In order to get the branching ratios, we compute the integral of the square of the fully 
symmetric amplitudes over the Dalitz plot (taking care to avoid sextuple counting). 
Doing this gives 

B(B + -> K + n + n-) fully sym = 0.19 B(B + K + n+n-) , 
B{B° d -> K + K°K-) fully sym = 0.50 B(B° d -+ K + K°K~) . (19) 

From Ref. [TO], we have r /r+ = 0.93, B(B+ ->■ K + n+n-) = (51.0 ±2.9) x lO" 6 and 
B(B° d -> K + K°K~) = (24.7 ± 2.3) x 10~ 6 . Eq. (US} then gives 

fsu(S) = 0-17 ±0.06. (20) 

(This error does not include the errors in the parameters obtained from the Dalitz- 
plot analyses of the two decays.) 

We can now put all the pieces together to describe how the fit is to be performed. 
The fully symmetric amplitudes for the decays B d — y K + -k q -k~ , B d — y K 7r + 7T~, 
B° d -> K + K°K~ and B° d ->■ K°K°K° are given in Eq. flU]). They are a function 
of 10 unknown parameters, including 7. The value of fsu(3) is taken from Eq. (120]) . 
The 11 observables and their errors are computed as described above - the (fully 
symmetric) branching ratios and direct CP asymmetries are given in Eq. (|T5|) . and 
the indirect CP asymmetries in Eq. ( TT6|) . Note that these are for specific values of 
s ++ and s + . One has a different set of observables for each (s ++ , s + ) pair. With 10 
unknowns and 11 constraints, one can now perform the fit. This will determine the 
magnitudes and relative strong phases of the five effective diagrams, as well as 7, 
all for the chosen values of (s++, s+). This is to be repeated for each independent 
(s ++ , s + ) paiid. This has two effects. First, one will be able to fix the momentum 
dependence of the diagrams. Second, and more importantly, since 7 is momentum 
independent, one can average over all the (s ++ ,s + ) fits. This will reduce its error, 
perhaps considerably. 

Now, we already have experimental information about most of the required B — > 
Kirn and B — y KKK decays. In particular, the measurements of the Dalitz plots 
of B° d -»■ K + n°n-, 5° -»■ K n + n- and 5° -»■ K + K K~ are described in Refs. [IT], 
|12j . and [H], respectively. On the other hand, we do not yet have the Dalitz plot 
of B d — y K°K°K°. The branching ratio and CP asymmetries of B d — y KsK$Ks 
are given in Ref. [17]. While the use of the final state KgKgKs is excellent - 

The two pairs (s++,s+)i and (s++,s+)2 are considered as independent if 
\M fully s V m((s++, s+)i)| and \M fully sym((s ++ ,s + )2)\ do not overlap when one takes into 
account the errors on the Dalitz-plot parameters of Eq. (TTU)) . 
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it is proportional to the fully symmetric state of K K K - the observables are 
momentum independent. That is, an integration over the Dalitz plot has been 
performed. However, the method described in this paper requires the momentum- 
dependent observables. Once the Dalitz plot for B® — > KgKgKs is known, this 
method for extracting 7 can be carried out. 

Even though all the experimental data is not yet available, we can still attempt to 
estimate the precision with which 7 can be obtained. Consider first B® — > K + K°K~ . 
According to the BaBar measurement in Ref . [H] , the largest contributions to this de- 
cay come from the <j)K° and foK° resonances, and the (K + K~) nrK°, (K + K°) nrK~ 
and (K~ K°)nrK + non-resonant pieces. They find 

c . = 0.0085 ± 0.0010 , 
c 3 = 0.622 ± 0.046 , 
Cj = 1 (fixed) , 
cj = 0.33 ± 0.07 , 

Cj = 0.31 ± 0.08 , (21) 

where Cj is defined in Eq. (flOj) . The errors, which are statistical only, range from 7% 
to 25%. The above method describes how to obtain Aif u ii y sym(B® — > K + K°K~) 
from the amplitude given in Ref. [9], and from this the B® — > K + K°K~ observables. 
A full numerical analysis is needed to do this, properly taking into account the errors 
on the Cj above, as well as the errors on the 8j and Fj of Eq. ( TTOj) . and the other 
resonances. However, a rough guess is that the errors on the observables will be 
about 20%. Similarly, we (guess)timate that the errors on the observables of the 
other decays, including those of B® — > K°K°K°, will be ~ 20%. In order to obtain 7, 
a fit to the observables must be performed, taking into account the SU(3)-breaking 
factor of Eq. ( 120]) (the error on fsu(3) will increase once the errors in the Dalitz- 
plot parameters are included), and one must average over the independent (s++, s + ) 
pairs. It is impossible to predict with any accuracy what the error on 7 will be, but 
an error of 0(25%) does not seem unreasonable. 

How does this compare with the precision on 7 measured in two-body decays? 
The answer is: not that badly. The standard way of directly measuring 7 uses 
B D°/D°K decays within the GLW [H] or ADS rj5] methods. The latest 
measurement yields 7 = (681^)° [IB] , i.e. the error is ~ 15%. To be sure, our 
estimated error of 0(25%) on the value of 7 as extracted from three-body decays 
is worse than 15%. However, it is still roughly the same size, and if a full analysis 
were done, the real error might turn out to be smaller than our estimate. More 
to the point, when the Dalitz-plot measurements are done at the super- 1? factory, 
the Dalitz-plot parameters will be obtained with a smaller statistical error. This 
will have two effects. First, the error on 7 will be reduced for each (s ++ , s + ) pair. 
Second, one will have more independent (s ++ , s + ) pairs, so the error will be further 
reduced when one averages over all the (s++, s + ) fits. (See the discussion following 



foK° 
(K + K-) NR K° 
(K + K°) NR K- 

(K~K°) \mK + 
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Eq. ( 1201) .) Thus, the extraction of 7 from three-body B decays may turn out to be 
more precise than that from two-body decays. 

Note added: after this paper was submitted, the Dalitz-plot analysis of B® — > 
KsKsKs was submitted to the arXiv, see Ref. [T7] . 
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